Real-space Renormalization Group Approach for the Corner Hamiltonian 
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We present a real-space renormalization group approach for the corner Hamiltonian, which is 
relevant to the reduced density matrix in the density matrix renormalization group. A set of 
self-consistent equations that the renormalized Hamiltonian should satisfy in the thermodynamic 
limit is also derived from the fixed point of the recursion relation for the corner Hamiltonian. We 
demonstrate the renormalization group algorithm for the 5 = 1/2 XXZ spin chain and show that 
the results are consistent with the exact solution. We further examine the renormalization group 
for the S = 1 Heisenberg spin chain and then discuss the nature of the eigenvalue spectrum of the 
corner Hamiltonian for the nonintegrable model. 



PACS numbers: 75.10.Jm, 05.50.+q 

I. INTRODUCTION 

Since Wilson's pioneering work for the Kondo im- 
purity problem, 0, |^ the numerical renormalization 
group (NRG) procedure has been developed extensively 
for the various low-dimensional quantum many-body sys- 
tems. In particular, the density matrix renormalization 
group(DMRG) enables us to obtain reliable results for 
the ground states of the one-dimensional(lD) quantum 
many-body systems. 0,0,1^ However, the DMRG has an 
essential difference from Wilson's approach; In the Wil- 
son's NRG, the lower-energy states of the Hamiltonian 
are selected to become relevant in low-energy physics, 
while the larger eigenvalue states of the reduced density 
matrix of the ground state are more important in the 
DMRG. Although the variational principle for the ma- 
trix product form of the ground state wavefunction lies 
behind the DMRG|^, what kind of the low energy theory 
is generated by the DMRG is not discussed from the Wil- 
son's renormalization group view. Of course, it may be 
difficult to answer such a fundamental question for gen- 
eral cases. In this paper, we want to present a possible 
approach to address the question. 

In order to discuss the reduced density matrix in the 
DMRG, it is essential to bring Baxter's corner transfer 
matrix(CTM), which is defined as the partition function 
of a quadrant for the 2D classical lattice systems. 0,11,13 
The peculiar feature of the CTM for the integrable model 
is that the CTM A can be exponentiated as 



A(/i) cx e 



(1) 



where /i is a spectral parameter and K is the corner 
Hamiltonian(the detailed definition is given by ©).[?', 
H, 1^ As pointed out in ref . jllj , the reduced density 
matrix in the DMRG is related to a product form of the 
CTMs. Then the reduced density matrix p for the in- 
tegrable model in the thermodynamic limit can be also 
written as 

P e-"^, (2) 

where a is a certain parameter. On the basis of eq. 

we can view that the DMRG provides a low-energy ef- 



fective theory of the corner Hamiltonian, suggesting that 
the real-space renormalization group approach also works 
for the corner Hamiltonian. Although eqs. iQJ or are 
established only for the integrable models, we may ex- 
tract implications about the role of the reduced density 
matrix for general cases from the corner Hamiltonian. 

In addition, the reduced density matrix in the DMRG 
is recently attracting much attention in the context of 
quantum information theory. 0, 0, 0| When a re- 
duced density matrix p is defined for a subsystem with 
a cut of certain geometry, the entanglement entropy 
S = — tr(/9lnp) plays an important role in understand- 
ing how quantum states are entangled through the other 
subsystem. The reduced density matrix or CTM in the 
DMRG is interestingly associated with the entanglement 
entropy of the ID quantum systems with the half-infinite 
cut, implying that the property of the corner Hamilto- 
nian is theoretically important from the quantum in- 
formation theory point of view. For instance, Peschel 
has studied a generalization of the corner Hamiltonian 
for the various free fermion models, which provides new 
insight in the studies of the properties of the quantum 
states. XL >1S, tlH 

In this paper, we discuss the real-space renormaliza- 
tion group approach for the corner Hamiltonian. In the 
next section, we formulate the recursion relation and the 
iterative diagonalization algorithm for the corner Hamil- 
tonian. We then derive the self-consistent equations for 
the renormalized corner Hamiltonian in the thermody- 
namic limit. We also rewrite the self-consistent equation 
so as to be more relevant to the block Hamiltonian in 
the DMRG formulation. In §3, we demonstrate the cor- 
ner Hamiltonian NRG for the XXZ spin chain. We find 
that the resulting low-energy eigenvalue spectrum of the 
corner Hamiltonian reproduces well the exact solution 
for the massive case. For the XXZ chain in the mass- 
less region, we analyze the size dependence of the NRG 
spectrum and verify the logarithmic size dependence ex- 
pected by the conformal field theory (CFT) [Hill. In §4, 
we further apply the NRG to the S = 1 Heisenberg spin 
chain, which is a typical example of the nonintegrable 
spin model. We then find that the dominant structure of 



2 



the eigenvalue spectrum of the corner Hamiltonian can 
be scaled by one parameter to the corresponding spec- 
trum(logarithm) of the reduced density matrix. More- 
over, we analyze the "spectral flow" of the corner Hamil- 
tonian for the S — 1 bilinear-biquadratic chain. In §5, we 
summarize our results and discuss further implications of 
the eigenvalue spectrum of the corner Hamiltonian. 



II. FORMULATION 

A. definitions 

In this section we consider the 5 = 1/2 XXZ spin chain 
for convenience. However, the arguments in the following 
can be generalized straightforwardly to the general ID 
quantum systems. We write the local Hamiltonian of the 
XXZ spin chain as 

where S is the 5=1/2 spin(not Pauli) matrix. The 
matrix element of Q is labeled by the spin indices s„, 
Sn+i and s^, s'^j^i- However, we do not show explicitly 
such spin indices without necessity. 

Let us denote the Hamiltonian of N spins as 

= ^ hn,n+l, (4) 

for which the free boundary condition is basically as- 
sumed. In the context of the DMRG, is the right- 
half(or left-half) block of the total Hamiltonian and, in 
this labeling of the site index, n = 1(N) is assigned to 
the center(edge) of half of the system. 

We further define the corner Hamiltonian as 

N-l 

-fS^N = ^ nhn,n+l, (5) 
n=l 

whose graphical representation is depicted in Fig. ^ As 
noted in the introduction, the corner Hamiltonian is the 
generator of the CTM for the integrable models. Thus, 
(O contains N{N — l)/2 local bonds, corresponding to 
the quadrant of the 6-vertex model. Here, we note that 
Kjsi is clearly an Hermite matrix. 



B. recursion relation 

In order to formulate the real-space renormalization 
group for the corner Hamiltonian, we have to set up the 
recursion relation for matrices having different dimension 
sizes. For this purpose, we introduce some notations for 
a 2^ X 2^ matrix X: 



O—O N 




1 2 3 N-l N 



FIG. 1: Graphical representation of corner Hamiltonian A"n. 
White circles indicate spins and the lines connecting two spins 
mean the local interaction /i„,„+i. The horizontal lines are 
"stacked" to become the corner Hamiltonian. 

where the row index of X-^ is labeled S2 • • • Sn-i-i and thus 
X^ has the index si,S2 ■ ■ ■ sn-i-i- Similarly, we also use 

^N-l = <5(Sl,s'l)^(s2,S2)^N-l• 
Wc then construct the recursion relation of the corner 
Hamiltonians between N and N-fl. As is illustrated in 
FigEl we can decompose Xn+i into three pieces and then 
find 

i^N+l = /ll,2 + i?N + ^^N- (7) 

However, this relation contains both of H and K, which 
is not convenient for capturing the eigenvalue structure 
of the corner Hamiltonian directly. In order to eliminate 
in lO, we exploit a supplemental recursion relation 
for which is represented as 

i/^ = XS-/fr_i. (8) 

We can thus construct the recursion relation consisting 
of the corner Hamiltonians 

i^N+i = /m.2 + 2XS-XS*_i. (9) 

A key point on ij^ is that we have derived the recursion 
relation between the corner Hamiltonians for three sizes 
(N-|-1,N,N-1) rather than for the two sizes (N-|-1,N) by 
eliminating H^. Here, it should be noted that such a 
construction of the recursion relation is almost parallel 
to the "logarithm" of Baxter's recursion relation for the 
CTMs0. 




1 2 — N+1 12 2 N+1 



"N+1 "1,2 

FIG. 2: Graphical representation of recursion relation (Q. 

We next convert the bases of the matrices into the 
representation diagonalizing K^: 



X^ = 5(si, s'O^N, 



(6) 



(10) 
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where cjn is a diagonal matrix of the eigenvalues of JCn, 
and J7n is a unitary matrix consisting of the correspond- 
ing eigenvectors. We assume that the elements of cjn are 
arranged in increasing order. Then, eq. ^ becomes 

i^N+i = h,,2 + 2C/MC^N - C^N-i^N-iC^N-r (11) 
Here, we introduce new unitary matrices 

P^ = U;*\U;, and Pn = C/;^C/n+i- (12) 

By using these matrices, we finally obtain the main re- 
cursion relation for the corner Hamiltonians 

i^N+l = i'N^/il,2PN + 2c^N - ^N^^N-l^N, (13) 

where 

^N+l = J/N+ii^N+lC^N+l = ^'n+I^^N+I^'n+i- (14) 

The above expression is very useful for the NRG calcula- 
tion, since the unitary matrix appearing in the right-hand 
side of eq. ((T^ is only P^. 



C. iterative diagonalization 

We solve the eigenvalue problem of the corner Hamil- 
tonian for a sufficiently large N recursively. We first set 
up the initial corner Hamiltonian for N=2 and 3 and 
diagonalize them. We then obtain an extended corner 
Hamiltonian by the recursion relation IjlSI) . and diago- 
nalize A'n+1 numerically by the Householder method, 

^n+i-Pn+1 = (15) 

By using these extended matrices, we can return to H13fl 
with incrementing N^N+1. 

In the above recursive step, the dimension of A'n+i 
is doubled. Thus, we retain half of the lower-energy 
eigenvalues, in each recursion step, and the correspond- 
ing eigenvectors Pn-i-i, after a sufficiently large N. If the 
number of retained eigenvalues is m, the dimension of 
^N+i is 2m. We then keep the lower-half eigenvalues, 
and Pn-i-1 becomes a 2m x m matrix, which plays the 
role of the projection matrix. In the context of retain- 
ing the lower-energy eigenvalues, the present algorithm 
is faithful to Wilson's original NRG method rather than 
the DMRG that keeps larger eigenvalues of the reduced 
density matrix . 



D. self-consistent equations 

After a sufficient number of iterations, the matrices will 
converge to those in the thermodynamic (N — > c») limit. 
In this sense, we omit the subscript N when discussing 



the matrices in N — > oo. Then we can write down a set 
of the self-consistent equations, 

K = P*^hi^2P* + 20;* - P*^uj**P*, (16) 
KP = Plu, (17) 
P^P = /, (18) 

which determine the eigenvalue structure of the corner 
Hamiltonian in the thermodynamic limit. Here, we note 
that the dimensions of the matrices are as follows: 

Lu : m X m, lu* : 2m x 2m, lo** : 4m x 4m 
P : 2m X m, P* : 4m x 2m, K : 2m x 2m. 

Unfortunately, we cannot illustrate an impressive graphi- 
cal expression for the above self-consistent equations, un- 
like those of the CTMs(see chap. 13 in ref. |^). However, 
we think that it is important to write down the closed 
form of the self-consistent equations in considering the 
fixed point of the DMRG. 

Further, we deduce self-consistent equations for the 
renormalizcd Hamiltonian, which may be relevant to the 
DMRG; Once the renormalizcd corner Hamiltonian is ob- 
tained, the renormalizcd Hamiltonian is also reproduced 
from © as 

H* =Lu* ~ P*''lu**P*, (19) 

where H* = U*'<H*U*. In connection with the DMRG, 
an alternative expression 

H ^ P*^hi^2P* - P*^uj**P*, (20) 

IS more suitable, where H = U*^ HU* . Substituting ^ 
into (fTB|l . we obtain 

K^H + uj*. (21) 

Further substituting this equation of K into (|17|l . we fi- 
nally arrive at 

HP = Puj- uj*P, (22) 

or equivalently, 

Y,H{s,,s[)P{s[) = [P{s,),cu]. (23) 

In the last equation, we recover the indices of the bare 
spins and thus (|23|l is the equation for the m x m matrices. 

Equations H18|) . 12U|I and H23|) . as a set, are also the 
equivalent self-consistent equations for the renormalizcd 
Hamiltonian in the thermodynamic limit. Here, it should 
be recalled that uj is the eigenvalue matrix not of the 
reduced density matrix but of the corner Hamiltonian. 
For the case of the integrable model, however, we have 
UJ (xlnpd with the common eigenvectors, where pd is the 
eigenvalue matrix of the reduced density matrix in the 
DMRG with an appropriate normalization. 

Of course the above-obtained self-consistent equations 
are too formal for practical use. In order to solve them, 
we have to employ the iterative method described in the 
previous subsection. 
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III. RESULTS FOR THE 5=1/2 XXZ CHAIN 

In this section, we demonstrate the NRG of the cor- 
ner Hamihonian for the S = 1/2 XXZ chain. We first 
discuss the Ising-like anisotropic case(A > 1), for which 
the corner Hamiltonian has a discrete eigenvalue struc- 
ture in the thermodynamic hmit. The NRG calculation is 
performed with the retained number of bases m = 200. 
Figure |3| shows the N dependence of the lowest seven 
eigenvalues for A = 2. In the figures, the "ground state" 
energy is set to be zero. For a small N, we can see the 
oscillation with respect to N=even or odd, for which the 
total iS*^ of the system takes an integer or a half integer 
alternately. As N increases, we can see that the eigenval- 
ues rapidly converge to the constant values corresponding 
to the thermodynamic limit. 




matrix. The coefficient a can be identified as 

4A 



sinh A ' 



(26) 



according to Baxter's parametrization of the Boltzmann 
weight of the 6-vertex model. Here, it should be noted 
that the above exact spectrum (|25|l with II24II was cer- 
tainly reproduced by the DMRG, as in Refs.|l2L[l^. 

Figure0]shows the lowest 100 eigenvalues of the corner 
Hamiltonian for the XXZ model with A = 2, which are 
obtained after a sufficient number of iterations. The cor- 
respondence between the results of the exact solution, the 
corner Hamiltonian NRG and the DMRG(2Qi] is perfect 
within the level of numerical accuracy. This implies that 
the NRG algorithm for the corner Hamiltonian works 
successfully. 




100 




100 



FIG. 3: Size dependence of lowest 7 eigenvalues of the corner 
Hamiltonian for the XXZ model with A = 2. The NRG 
calculation is performed with m — 200. The ground state(i — 
1) corresponds to the horizontal axis. The i = 3, 5 and i = 5, 6 
states are two-fold degenerating. 

The exact eigenvalue structures of the corner Hamil- 
tonian and the reduced density matrix of the XXZ 
model in the thermodynamic limit are well known for 
A > 1, where the parameterization A = cosh A is very 
useful. 0, II, la| Through the diagonal representation of the 
GTM, we can extract the exact spectrum of the corner 
Hamiltonian as 

u; = . diag(0, 0, 1, 1, 2, 2, 3, 3, 3, 3 • • •), (24) 

where diag(- • •) denotes a diagonal matrix whose diagonal 
entry is • • •. Since the wavefunction of the XXZ chain 
is represented as ^(A) oc A{X — fi)A{fi) on the basis of 
the corresponding 6-vertex model, we also have the exact 
spectrum of the reduced density matrix 

Pd{X) cx e-"", (25) 

where pd is the eigenvalue matrix of the reduced density 



FIG. 4: Lowest 100 eigenvalues of the corner Hamiltonian 
for the XXZ model with A = 2 for N = 200. The NRG 
calculation is performed with m = 200. 

Next we discuss the XXZ model with A < I, which has 
the gapless ground state. It is well known that the corre- 
sponding GTM is not normalizable in the thermodynamic 
limit. In other words, — » as N increases. However, 
the size dependence of the spectrum of the GTM or the 
corner Hamiltonian can be analyzed in the framework of 
the finite size scaling with the CFT; the conformal map- 
ping from the upper half-plane into the helical stairway 
geometry yields the spectrum of the GTM with an in- 
finitesimal transfer angle, which can be associated with 
the corner Hamiltonian. |2lL 122 . l2a | However, actual nu- 
merical calculations of the corner Hamiltonian for a large 
but finite size system have been difficult, away from the 
free fermion model. 

We define the energy gap of the corner Hamiltonian as 
Aw = (jj(^) — cj'^"^ for N =even, where the superscripts 
with the bracket denote a "quantum number" for the 
energy level taking account of the degeneracy structure; 
e.g., Au! corresponds to the energy difference between the 
first step and ground step of the stairwaylike spectrum. 
Then, the GFT argument yields the size dependence of 
Aw as 



Alu cx 



1 



ln(N/c) ' 



(27) 
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where c is an ultraviolet cutoff parameter. 



IV. APPLICATION TO S = 1 CHAINS 




FIG. 5: Size dependence of gap Auj for the XXZ cliain in tfie 
gapless regime up to N~4000. Tlie solid symbols indicate the 
NRG results. The lines are the fitted curves with eq. 1281 . 
The data correspond to A = 1.0, 0.5 and 0.0 from top to 
bottom. 

In Fig. 13 we show the size dependences of Aui for A = 
0.0, 0.5 and 1.0, which are obtained with m = 1000. We 
have confirmed that the m dependence of Auj is negligible 
up to N '--^ 4000. In order to check H27(l . we examine a 
fitting of the form: 



Ci 



InN + Ca' 



(28) 



where Ci and C2 are fitting parameters. Here we note 
that the parameter C2 is definitely requested, since a 
variable in a logarithm should be dimensionless. In Fig. 
El we can see that the NRG data is well fitted by We 
show the obtained values of Ci in Fig. The coefficient 
Ci may be connected to the scaling dimension in CFT. 
For this purpose, however, we need to carefully calculate 
the velocity with another numerical method. A further 
study is clearly desired for the analysis of the coefficient 
Ci. 



Let us examine the corner Hamiltonian NRG for the 
5 = 1 Heisenberg spin chain, which is a typical exam- 
ple of nonintegrable models. Since the S = 1 Heisenberg 
chain is in the gapful phase, the corner Hamiltonian and 
the reduced density matrix exhibit discrete structures in 
their eigenvalue spectra. Indeed we have confirmed that, 
as N increases, the low-energy eigenvalues rapidly con- 
verge to those in the thermodynamic limit. 

In order to compare them systematically, we define the 
logarithm of the eigenvalue spectrum of the reduced den- 
sity matrix as 



a 



(29) 



where we assume that the largest eigenvalue of pd is nor- 
malized to be unity. Here, we should recall that the direct 
relationship between the reduced density matrix and the 
corner Hamiltonian such as eq. ||2J) is not established for 
the non-integrable model, implying that the coefficient a 
is unknown a priori. Thus, we set a = — In p'"J^ / u^^'^ ; we 
scale WD so that its excitation gap corresponds to that of 
w, and then compare their stairwaylike structures. 



S=l Heisenberg chain 



o CO 
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100 




FIG. 6; Fitting resuh of Ci in ^ for the XXZ chain of 
< A < 1. The line is a guide for eyes. 



FIG. 7: Eigenvalue structure of the corner Hamiltonian for 
the S = 1 Heisenberg model. 

In Fig. [7| we show the comparison between the NRG 
result for the corner Hamiltonian and the coo spec- 
trum obtained with the DMRG, for which we have used 
a = 0.3942. Both the NRG and DMRG calculations are 
performed with m = 200. Since the integrability does 
not hold for the 5* = 1 Heisenberg chain, we can see that 
the regular stairway structure is disturbed. However, the 
remarkable point is that the dominant structures of the 
spectra for the corner Hamiltonian and the reduced den- 
sity matrix are very similar; Since is scaled to to by 
one parameter a, the NRG for the corner Hamiltonian 
can contain almost the same level of information as the 
DMRG for the S = 1 Heisenberg model. In ref.^^, we 
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have found that the eigenvalue structure of the reduced 
density matrix seems to have universal asymptotic behav- 
ior not only for the integrable models but also for a class 
of ID non-integrable models. The above one-parameter 
scalability with respect to a may be related with such 
universal asymptotics of the eigenvalues. 

In order to clarify the nature of the corner Hamiltonian 
further, we investigate the 5 = 1 bilinear-biquadratic 
chain 



(30) 



where S is the S = 1 spin matrix. This model has 
been intensively studied in the context of the Haldane 
gap. /3 = ±1 is exactly solvable and the ground state 
is gapless.m /3 = 1/3 is the AKLT model whose 
ground state is exactly represented as the valence-bond- 
solid(VBS) state, which is a typical example of the matrix 
product form of the wavefunction. [2^ From the DMRG 
point of view, the AKLT model is particularly important, 
since the DMRG with m = 4 can exactly reproduce the 
VBS state. Thus, it is interesting to know how the 
spectrum of the corner Hamiltonian for (|30|l behaves, in 
contrast with the DMRG. 

In Figure |S1 we show the /3-dependence of the lowest 
64 eigenvalues for -0.8 < (3 < 0.8. For -0.5 < (3 < 0.8, 
the NRG iteration sufficiently converges to the thermo- 
dynamic limit within m = 400 and N=400. As /3 ap- 
proaches ±1, the convergency becomes worse, since the 
exactly solvable points(/9 = ±1) are gapless. Thus, we 
improve the accuracy up to m = 800 and N=2000 for 
(3 < —0.6. In the vicinity of the gapless points(/3 < —0.8 
or /3 > 0.8), unfortunately, we do not succeed in extract- 
ing the spectnmr of the thermodynamic limit. 
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FIG. 8: /3-dependence of lowest 64 eigenvalues of the corner 
Hamiltonian for the S = 1 bilinear-biquadratic model. An 
open circle corresponds to a step of the almost degenerat- 
ing eigenvalues in the stairwaylike eigenvalue structure. The 
arrow indicates the AKLT point(/3 = 1/3). 

In Fig. |S1 we see only a few points of levels for a fixed 
/3, since the corner Hamiltonian has the almost degen- 
erating eigenvalue structure illustrating the stairwaylike 
spectrum as in Fig. \7\ An interesting point is that the 



rearrangement of the spectrum clearly occurs, as (3 varies 
from —0.8 to 0.8. For instance, we can see that the 2nd 
and 3rd excitation levels merge as ^ 1. In the vicinity 
of P — ±1, the eigenvalue spectrum is expected to have 
a nearly regular structure governed by the integrability 
of the massless fixed point and the most relevant per- 
turbation associated with f3, although the corner Hamil- 
tonian at the massless fixed point itself is not normal- 
izable. The parameter /3 connects these two integrable 
limits continuously. Thus, the dominant structure of the 
spectrum is reconstructed adiabatically as /3 varies be- 
tween — 1</3<1. The present result demonstrates 
such reconstruction of the structure of the spectrum. 

Another important aspect of Fig. |S1 is that /? = 1/3 is 
not any special point for the corner Hamiltonian. Since 
the two-body Hamiltonian H30() at /? = 1/3 is the pro- 
jection operator for the S = 2 state on the nearest- 
neighboring bond, the "ground state" of the corner 
Hamiltonian itself is identical to the VBS state for the 
AKLT Hamiltonian. Accordingly, the convergence of the 
NRG iteration for (3 = 1/3 is very fast; after a few it- 
erations, the eigenvalue spectrum reaches the thermo- 
dynamic limit. However, we find that the spectrum 
of the corner Hamiltonian behaves continuously around 
(3 = 1/3; the dimension of the corner Hamiltonian clearly 
increases with respect to the system size, while the eigen- 
values of the reduced density matrix for the VBS state are 
exactly truncated by m = 4. This implies that the corner 
Hamiltonian certainly involves information on the exci- 
tations(in the sense of the usual Hamiltonian), while the 
reduced density matrix is determined only by the ground- 
state wave function. A further analysis of the deviation 
between these two spectra around /? = 1/3 may provide 
important information for understanding the correlation 
effects in the ID quantum systems. 



V. DISCUSSION 

In this paper, we have discussed the real-space renor- 
malization group for the corner Hamiltonian. We have 
derived the recursion relation for the corner Hamiltonian 
and its self-consistent equations in the thermodynamic 
limit, where the lower-energy states of the corner Hamil- 
tonian are more relevant. For the integrable model, we 
have verified that the corner Hamiltonian NRG yields 
the equivalent of the exact eigenvalue spectrum of the re- 
duced density matrix. Indeed we have shown that the ex- 
act spectrum is reproduced for the XXZ chain of A > 1. 
For the massless case, the size dependence of the excita- 
tion gap of the corner Hamiltonian is consistent with the 
logarithmic behavior expected by the CFT, although a 
further study is desired for the direct connection of the 
fitting results with the scaling dimension. For the 5 = 1 
Heisenberg chain, we have found that the dominant struc- 
ture of can be scaled to lo by one parameter a. More- 
over, the spectral flow for the bilinear-biquadratic chain 
exhibits a clear reconstruction of the spectrum between 
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the two integrable points (/3 = ±1) 

As far as the relation ^ is estabUshed, the present 
results imply that we do not need the diagonalization of 
the super block Hamiltonian(in the DMRG sense). This 
suggests that a further simplification of the DMRG algo- 
rithm might be possible in principle. Of course, we still 
have a gap to the practical use of the corner Hamiltonian 
NRG for calculating the physical quantities of the non- 
integrable cases, where the parameter a is not known a 
priori. However, it is important that the corner Hamilto- 
nian involves the same level of implications as the DMRG 
even for the non-integrable case. 

From the theoretical viewpoint, we have seen that the 
renormalization group for the corner Hamiltonian pro- 
vides a mount of interesting problems. It is intriguing to 
know how we can use the self-consistent equations 
((T7jl . lfTS|l or lO, (1201), (EH for an analysis of the ma- 
trices in the thermodynamic limit. For the critical sys- 
tem, the finite size analysis of the spectrum is explained 
by the CFT; the direct estimation of the scaling dimen- 
sion from the fitting data is a remaining problem. For 
the nonintegrable case, we have seen the two aspects of 
the corner Hamiltonian: the one-parameter scalability 
for the dominant structure of the spectra between the 
corner Hamiltonian and the reduced density matrix for 
the S* = 1 Heisenberg chain and, in contrast, their devi- 
ation for the AKLT model(/3 = 1/3). In addition to the 
above, we should remark on the boundary condition of 
the corner Hamiltonian; Clearly the corner Hamiltonian 



treated here is only the right(or left) block in the DMRG 
and the free boundary condition is imposed both for the 
center and edge spins. In contrast, the center spins in 
the DMRG are directly entangled with the other block 
of the system. This may be a key point in discussing the 
role of the entanglement of the states in the ID quan- 
tum systems, where a may have an important physical 
meaning even for the nonintegrable case. We also note 
that the connection to the density-matrix spectrum in 
the finite size DMRG, which we have not treated here, is 
a remaining important problem from the practical point 
of view. 

Our motivation at the start is to answer the question: 
what kind of low-energy effective theory is obtained by 
the DMRG in the Wilson's sense. We think that the 
present corner Hamiltonian approach provides a possi- 
ble way to address the question and, at the same time, 
stimulates further investigations on the issue. 
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